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Abstract 

Homogenized composite materials (HCMs) can support a singular form of optical propagation, known 
as Voigt wave propagation, while their component materials do not. This phenomenon was investigated 
for biaxial HCMs arising from nondissipative isotropic dielectric component materials. The biaxiality of 
these HCMs stems from the oriented spheroidal shapes of the particles which make up the component 
materials. An extended version of the Bruggeman homogenization formalism was used to investigate the 
influence of component particle orientation, shape and size, as well as volume fraction of the component 
materials, upon Voigt wave propagation. Our numerical studies revealed that the directions in which 
Voigt waves propagate is highly sensitive to the orientation of the component particles and to the volume 
fraction of the component materials, but less sensitive to the shape of the component particles and less 
sensitive still to the size of the component particles. Furthermore, whether or not such an HCM supports 
Voigt wave propagation at all is critically dependent upon the size of the component particles and, in 
certain cases, upon the volume fraction of the component materials. 
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1 Introduction 

A composite material, comprising random distributions of two (or more) particulate component materials, 
may be regarded as effectively homogeneous provided that the particles which constitute the component 
materials are much smaller than the wavelengths involved [1, 2]. Depending upon the shape, size and 
distribution of the component particles, such homogenized composite materials (HCMs) can exhibit optical 
properties not exhibited at all by their component materials, or at least not exhibited to the same extent by 
their components. Indeed, through judicious design, HCMs arising from commonplace component materials 
may be conceptualized which support rather exotic (and potentially-useful) optical phenomenons, while their 
component materials do not. A prime example is provided by HCMs which support plane- wave propagation 
with negative phase velocity [3, 4]. Another example — which provides the setting for the present study — is 
furnished by HCMs that support a singular form of optical propagation known as Voigt wave propagation 
[5]. 

In order to describe what constitutes a Voigt wave, it is helpful to first consider the nonsingular case of 
optical propagation in a linear, homogeneous, anisotropic, dielectric material. Generally two independent 
plane waves, with orthogonal polarizations and different phase velocities, can propagate in a given direction 
[6]. However, as first reported by Voigt in 1902 [7] and later described by others [8, 9, 10, 11, 12], in certain 
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dissipative biaxial crystals there are particular directions along which these two waves coalesce to form a 
single plane wave. A prominent feature of this coalescent Voigt wave is that its amplitude has a linear 
dependence upon propagation direction. Bianisotropic materials offer greater scope for Voigt waves [13], but 
herein our attention is restricted to the simpler case biaxial dielectric materials. 

Previously, the standard Bruggeman homogenization formalism was used to establish that certain dissi- 
pative biaxial HCMs can support Voigt wave propagation [5]. A follow-up study based on the second-order 
strong-permittivity- fluctuation theory — which represents a higher-order formulation of the standard Brugge- 
man formalism wherein two-point statistical correlations between particles of the component materials are 
taken into account [14] — emphasized the importance of correlation length for Voigt wave propagation 
[15]. Both of these earlier studies concerned HCMs arising from two uniaxial dielectric component materials 
(which cannot themselves support Voigt wave propagation). However, a biaxial dielectric HCM can also 
arise from isotropic dielectric components in instances where the shapes of the component material particles 
are non-spherical. For example, if each of the two component materials comprises an assembly of oriented 
spheroidal particles then the corresponding HCM will be biaxial, in general [16]. This is the scenario ex- 
plored here. We implement an extended version of the standard Bruggeman formalism [17], which utilizes a 
recently-developed extended depolarization dyadic formalism [18] in order to take into account the nonzero 
size of component particles. This approach allows us to investigate the influence of the non-electromagnetic 
attributes of the component materials — those being the orientation, shape and size of the component 
particles as well as the volume fraction of the component materials — upon the propagation of Voigt waves. 

In the notation adopted, vectors are represented in boldface, with the " symbol denoting a unit vector. 
Thus, the unit Cartesian vectors are written as x, y and z. Double underlining with normal typeface signifies 
a 3x3 dyadic; the identity 3x3 dyadic is / = xx + yy + zz; and the superscript T denotes the dyadic 
transpose. Double underlining with blackboard bold typeface signifies a 6x6 dyadic. The permittivity and 
permeability of free space are written as eo and /io, respectively. The free-space wavenumber is ko = uj-yje^jl^^ 
with uo being the angular frequency. 

2 Homogenization formalism 
2.1 Component materials 

We consider optical propagation in an HCM derived from two component materials, labelled a and h. The 
component materials are both taken to be isotropic dielectric materials with permittivity dyadics e = 606^/ 
and = Co 65/. The volume fraction of material a is fa while that of material b is fb = 1 — fa. 

Each component material comprises a randomly-distributed assembly of spheroidal particles. All material 
a particles have the same orientation and all material b particles have the same orientation, but these two 
orientations are generally different. For simplicity, both material a and b particles are assumed to have the 
same shape. The surfaces of the component particles, relative to their centres, are prescribed by the vector 



Here, r is the radial vector prescribing the surface of the unit sphere, the real-symmetric surface dyadic 
maps the spherical surface onto a spheroidal one, and 7^^ > is linear measure of particle size. In conformity 
with the homogenization regime, rji must be much smaller than the wavelengths involved, but — unlike in 
conventional approaches taken in homogenization studies [19] — we shall not insist that r]£ is vanishingly 
small. Furthermore, let us assume that r]a = Vb^ and henceforth simply write r] in lieu of r]£. 

The symmetry axis of the spheroidal particles comprising material a is taken to lie in the xy plane at an 
angle (p to the x axis. Thus, the surface dyadic for material a may be expressed as 



r£ = r]ilLp' ^ 



{£ = a,b). 



(1) 




+ /7(yy + zz)]-^^((^), (t/, t/. > 0) 



(2) 



R {(p) = cos (/?(xx + yy) +sin(/9(xy — yx) + zz. 



(3) 
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Without loss of generality, the material b particles are assumed to be aligned with the x axis. Thus, the 
surface dyadic for material b may be expressed as 



[[/,xx + [/(yy + zz)], (t/,t/.>0). 



(4) 



2.2 Homogenized composite material 



As the symmetry axes of the spheroidal particles comprising components a and b are not generally aligned, 
the resulting HCM is a biaxial dielectric material with a symmetric permittivity dyadic of the form 



-HCM 



: eo [e^xx + e^^yy + et (xy + yx) + e^zz]. 



(5) 



We implement an extended version of the Bruggeman homogenization formalism in order to estimate Ljj^j^- 
Accordingly, i^^^ is extracted from the nonlinear dyadic equation [20] 

/"{(i.-I^Cm)- [l + ^a-(ia-|^/CM)]''}+^{(l6-|^/CM)- [| + ^6-(|6-|^/Cm)]"'}=0' 

(6) 

using standard numerical techniques, such as the Jacobi method [21]. The depolarization dyadics ^ herein 
may be regarded as sums of two terms; that is. 



D=D°+D+ {£ = a,b). 



(7) 



The term D^^ represents the depolarization contribution arising from a vanishingly small particle described 
by the surface dyadic as given by the double integral [22, 23] 



-I pZTT pTT (J 

D° = — I / — ^ 



= a,b), 



(8) 



wherein the unit vector q = sin^cos(/)x + sin^sin(/)y + cos6>z. The depolarization contribution which 
derives from the nonzero size of the component particles is provided by the term D_'^ . It is most conveniently 
expressed in terms of elements of the 6x6 dyadic D"^ , per 



(m,nG {1,2,3}), 



(9) 



with [18] 

111)+ — 



47r/io 



sm( 



Je=o 



1 



K-^ — K 

+det 



{UTe' ■ q) ■ gHCM ■ ^7' ■ q)] ■ q) • fe7' • q) 



2q^ 



det rA(0)l ^ 



dO ^ 



{£ = a,b). 



Herein k,± are the q'^ roots of det \^{U_ ^ • q)] = 0, the vector q = ^q, while the 6x6 dyadics 



(p) 



-HCM 



(p/cj) X / 



(10) 



(11) 
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and 

£^(p)=A-\p)- lim A-^(p). (12) 

— IpI^oo — 

Analytical evaluations of the integrals in eqs. (8) and (10) are available for relatively simple anisotropic 
HCMs [24, 25], but for general biaxial HCMs numerical methods are needed to evaluate these integrals. 



3 Voigt wave propagation 



Let us turn now to the possibility of Voigt wave propagation in the HCM. All propagation directions relative 
to the symmetry axes of the HCM should be considered. It is expedient to do so indirectly, by investigating 
Voigt wave propagation along the z axis for all possible orientations of the HCM. Thus, we introduce the 
HCM permittivity dyadic in the rotated coordinate frame 



-HCM 



(a,/3,7) 



-HCM 

eiixx + e22yy + e33zz + ei2 (xy 
+ei3 (xz + zx) + e23 (yz + zy). 



c(^)-i;(/^)-c(7) 



-yxj 



where the orthogonal rotation dyadic 

R {P) = cos /3 ( X X + z z ) + sin /3(zx — xz) + yy. 



(13) 
(14) 

(15) 



and a, /3 and 7 are the three Euler angles [26]. 

In order for Voigt waves to propagate along the z axis of the biaxial dielectric material described by the 
permittivity dyadic (13), the following two conditions must be satisfied [27]: 

(i) V(a,/3,7) = 0, and 

(ii) W{a,f3n)^0. 
where the scalars 

V(a,/3,7) = 



13 ■ 



- 623 - 2623^33 [2612613 - (en - 622 ) 623] + (611-622) +46^2 



"^33 



and 



+2613 { ^23^13 - [2612623 + ( 611 - 622 ) 613 ] 633 } 



W(a, /3, 7) = 612633 - 613623 . 



(16) 



(17) 



Let us note that the conditions (i) and (ii) cannot be satisfied by isotropic or uniaxial dielectric materials. 



4 Numerical studies 
4.1 Preliminaries 

We now investigate the scope for Voigt wave propagation in the biaxial HCM, in terms of the shape, size and 
orientation of the component material particles, as well as the volume fraction of the component materials, 
by means of representative numerical calculations. Two stages are involved: first, 6^^^ is estimated using 
the extended Bruggeman formalism; and second, the quantities F(a, /3, 7) and W(a, ^, 7) are calculated as 
functions of the Euler angles. More specifically, the angular coordinates (a,/3,7) of the zeros of and the 
corresponding values of are computed. Notice that the angular coordinate 7 may be eliminated from our 
investigation because propagation parallel to the z axis (of the rotated coordinate system) is independent of 
rotation about that axis. 
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In the following, we chose = 1.5 and 65 = 12 as the relative permittivities of the component mate- 
rials a and h. Since ea,6 G R, the component materials are nondissipative. However, under the extended 
Bruggeman homogenization formalism with 77 > 0, the relative permittivity parameters of the HCM are 
complex- valued. The imaginary parts of the HCM's relative permittivity parameters are indicative of losses 
due to scattering from the macroscopic coherent field [28]. This is a general feature of higher-order ap- 
proaches to homogenization, as occurs in a similar fashion with the strong-permittivity-fluctuation theory 
[29], for example. 

We take the spheroid parameters Ux = ^ ^ p and U = 1 — (p/18), and consider the range < p < 9. 
Thus, the component particles become increasingly elongated as the eccentricity parameter p increases from 
zero, whereas in the limit p ^ the particle shape becomes spherical. 

4.2 HCM constitutive parameters 

The extended Bruggeman estimates of the HCM's relative permittivity parameters ex,y,z,t are plotted as 
functions of spheroid orientation angle (f and volume fraction fa in Fig. 1. Here, we fixed the relative size 
parameter r] = 0.2 /ko and the eccentricity parameter p = 9. We see that the real parts of Cx^y^z decay almost 
linearly as fa increases from to 1, but Ke{ex,y,z} are largely insensitive to variations in (f. While the real 
part of et is nonexistent for (p = 0° and 90°, and also in the limits fa^O and 1, away from these boundary 
values we have Re{et} < with a local minimum occurring at approximately (f = 45° and fa = 0.26. The 
imaginary parts of Cx^y^z attain their largest values in the vicinity of fa = 0.26, and vanish at = and 
1. As (p increases from zero, Imjea^} generally increases, Imje^y} generally decreases, and Imje^} is largely 
unchanged. The imaginary part of is null- valued along the boundaries cp = 0° and 90°, and /a = and 
1, but away from these boundary values Im{et} exhibits a local maximum which coincides with the local 
minimum exhibited by Re{et}. The following two limits, which hold for arbitrary fa G (0, 1), are especially 
noteworthy: (a) As ^ 0°, we find that Cx ^ Cy = Cz] i.e., the HCM is uniaxial. Accordingly the HCM 
cannot support Voigt wave propagation in this hmit [27]. (b) As (p ^ 90°, we find that Cx, Cy and are all 
different while ct = 0; i.e., the HCM's biaxial structure is orthorhombic. 

Next we turn to the dependency upon the component particle shape and size. In Fig. 2, ex,y,z,t are 
plotted against relative size parameter kof] and the eccentricity parameter p. Here, we fixed the volume 
fraction fa = 0.25 and the spheroid orientation angle cp = 45°. The real parts of ex,y,z generally increase 
as 77 increases, whereas Re{et} is largely independent of 77. As p increases, Rejea;} generally increases 
but Rele^y,^,^} generally decrease. The most conspicuous feature of the plots of lm{ex,y,z,t} is that these 
quantities increase uniformly from zero as rj increases from zero. Also, Im {cx^y^z} are fairly insensitive to 
increasing p but Im{et} increases markedly as p increases. Let us note that in the limit p ^ 0, we have 
= ^y = and Ct = 0; i.e., the HCM is an isotropic dielectric material, regardless of the size parameter r] 
(or the orientation angle (p and volume fraction fa). 

4.3 Orientations for Voigt waves 

Before presenting in detail the results of our numerical study of Voigt wave propagation in the HCM char- 
acterized in Figs. 1 and 2, the following should be pointed out. In general, for a given HCM the Voigt wave 
conditions Y = and W ^ are satisfied at two distinct orientations, specified by the angular coordinates 
<^ = <^i,2 and P = Pi ^2- However, for the numerical investigations reported herein we found that ai ^ a2 
and /^i ~ /32, with the differences between the two orientations being at most approximately 1% and often 
much less. Consequently, each curve of a, corresponding to a zero of IF], presented in the following Figs. 3, 
6 and 7 would appear at higher resolution as two very closely-spaced curves, namely that of ai and that of 
a2. And similarly for each curve of /3, and the corresponding curves of in Figs. 3, 6 and 7. However, at 
the resolution of Figs. 3, 6 and 7 this is difficult to perceive. This matter is elaborated upon in due course, 
in the discussion of Figs. 4 and 5. 

Let us now begin by considering the effects of volume fraction /«, and spheroid orientation angle cp^ on 
Voigt wave propagation in the HCM characterized in Figs. 1 and 2. The angular coordinates a and P of 
the zeros of along with the corresponding values of are plotted versus fa G (0,1) in Fig. 3 for 
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if = 30°, 60° and 90°. The a coordinate for if = 90° changes abruptly at fa = 0.48, from being very close 
to 90° for fa < 0.48 to being very close to 180° for fa > 0.48. The change in the a coordinate for mid-range 
values of fa becomes progressively less abrupt as cp decreases: for cp = 60° the corresponding graph of a 
is approximately sigmoidal while for cp = 30° the graph is nearly linear. The /3 coordinate is also highly 
sensitive to the volume fraction. Indeed, for cp = 90°, the value of P ranges from 0° at fa = 0.48 to 90° in 
the limits fa^O and 1. This sensitivity becomes less pronounced as (p decreases. From the corresponding 
graphs of we can deduce for which values of fa the HCM supports Voigt wave propagation. In the 
limits /a ^ and 1 we see that ^ and therefore the HCM cannot support Voigt waves in these limits. 
There are three further points for (p = 90°, namely fa = 0.2, 0.48 and 0.83, at which \W\ = and the HCM 
cannot support Voigt wave propagation. 

As mentioned earlier, each curve in Fig. 3 (and indeed in Figs. 6 and 7) actually represents two closely- 
spaced but distinct directions for Voigt wave propagation. In order to better appreciate this feature, in Fig. 4 
the curves of ai and a2 for cp = 90° in Fig. 3 are reproduced at much higher resolution. The corresponding 
curves of /3i^2 are indistinguishable even at this higher resolution. At values of fa < 0.48, the graphs of ai 
and a2 are mirror images with respect to reflection about ai^2 = 90°, and the same applies to the ai^2 curves 
for fa > 0.48 but with respect to the symmetry hue ai^2 = 180°. The differences between the ai and a2 
curves ranges from approximately 1° down to 0°. Notice that the fa values at which the difference is 0°, i.e., 
fa = 0.2 and 0.83, correspond to the zeros of |W| observed in Fig. 3. 

In order to illustrate the two distinct orientations for Voigt wave propagation for cp ^ 90°, in Fig. 5 
a representative example is provided for cp = 30°, with fa = 0.3 and the eccentricity and size parameters 
being the same as in Fig. 3. Here the normalized values of |y| are plotted versus the angular coordinates 
a e (164.5°, 165.5°) and f3 G (70.4°, 71.4°). Two \Y\ zeros can be seen: one at a = 164.8°, f3 = 70.9° and 
the other at a = 165.2°, f3 = 70.9°. 

Next the influence of the component particles' shape on the propagation of Voigt waves is considered. 
The angular coordinates a and f3 of the zeros of IF], along with the corresponding values of are 
plotted as functions of the eccentricity parameter p G (0,9) in Fig. 6 for the spheroid orientation angle 
(p G {30°, 60°, 90°}. The volume fraction was fixed at fa = 0.25 and the size parameter 77 = 0.2//co- The 
values of a for each value of cp are quite different, but these values are almost independent of p. The values of 
P for each value of cp are also quite different, but for the values of this angular coordinate decrease gradually 
as p increases. The corresponding values of increase as p increases, most rapidly for cp = 90° and least 
rapidly for cp = 30°. Furthermore, ^ as p ^ 0, regardless of cp^ as would be expected since the HCM 
becomes an isotropic dielectric material in this limit. 

Lastly, we turn to the effect of the component particle size, as accommodated by the extended Bruggeman 
formalism. In Fig. 7, graphs of the angular coordinates a and P at which \Y\ = 0, and the corresponding 
values of against the relative size parameter kof] G (0,0.2) are displayed for the spheroid orientation 
angle (p G {30°, 60°, 90°}. Here the volume fraction fa = 0.25 and the eccentricity parameter p = 9. While 
the size parameter r] has very little influence upon the orientations for Voigt waves, as indicated by the 
nearly horizontal graphs of a and /3, it does have a profound influence upon whether or not Voigt waves can 
propagate. Since \ W\ ^ as 7^ ^ 0, we infer that the nondissipative HCM — arising from vanishingly-small 
component particles — does not support Voigt wave propagation. In contrast, for nonzero values of r] we 
see that > and therefore Voigt wave propagation is supported. 

5 Closing remarks 

Engineered materials, in the form of HCMs, may be conceptualized which support Voigt wave propagation 
while their component materials do not. The case considered here involved remarkably simple component 
materials, namely nondissipative isotropic dielectric materials, in contrast to earlier studies on Voigt- wave- 
supporting HCMs which involved dissipative uniaxial dielectric components [5, 15]. In the present case, 
the ability of the HCMs to support Voigt waves relied upon the shape, orientation and nonzero size of the 
component particles, as well as the volume fraction of the component materials. In order to cater for such 
component materials, an extended version [18] of the well-established Bruggeman homogenization formalism 
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[17, 19] was needed. We note that the homogenization approaches adopted in ear her Voigt wave studies, to 
wit the Maxweh-Garnett formahsm [5] and strong-permittivity-fluctuation theory [15], could not be used here 
because they cannot accommodate component particles with differing orientations or component particles 
of nonzero size. 

Our numerical investigations revealed that the directions in which Voigt waves may propagate is highly 
sensitive to the orientation of the component particles and to the volume fraction of the component materials, 
but less sensitive to the shape of the component particles and less sensitive still to the size of the component 
particles. Furthermore, whether or not a HCM supports Voigt wave propagation is critically dependent 
upon the size of the component particles and, in certain cases, upon the volume fraction. For the scenarios 
considered at here, it was found that, in general, the two directions which support Voigt wave propagation 
are very close together. This is consequence of the component materials being nondissipative. For dissipative 
component materials, these two directions can be more widely spaced [15]. 

This study further emphasizes the role of the micro- and/or nano-structure in determining the macroscopic 
optical properties of engineered materials, and paves the way for a study of Voigt wave propagation in 
bianisotropic HCMs, which — courtesy of their enormous parameter space — present greater scope for Voigt 
waves. 
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Figure 1: The extended Bruggeman estimates of relative permittivity parameters of the HCM plotted versus 
volume fraction fa G (0, 1) and spheroid orientation angle G (0°, 90°). The relative size parameter kof] = 0.2 
and the eccentricity parameter p = 9. 
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Figure 2: As Fig. 1 except that the HCM's relative permittivity parameters are plotted versus the relative 
size parameter koT] G (0,0.2) and the eccentricity parameter p G (0,9). The volume fraction fa = 0.25 and 
spheroid orientation angle cp = 45°. 
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Figure 3: The angular coordinates a, /3, and the absolute value of the quantity plotted versus volume 
fraction fa for spheroid orientation angle = 90° (blue, dashed curves), 60° (green, solid curves) and 30° 
(red, broken dashed curves). The size parameter r] = 0.2//co and the eccentricity parameter p = 9. 
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Figure 4: Close-up plots of a versus fa in Fig. 3 for = 90° 
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Figure 5: The normalized values of \Y\ plotted versus angular coordinates a G (164.5°, 165.5°) and j3 G 
(70.4°, 71.4°) for spheroid orientation angle = 30°, volume fraction fa = 0.5, eccentricity parameter p = 9 
and size parameter r] = 0.2/ko. 
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Figure 6: The angular coordinates a, /3, and the absolute value of the quantity VK, plotted versus the 
eccentricity parameter p for spheroid orientation angle (p = 90° (blue, dashed curves), 60° (green, solid 
curves) and 30° (red, broken dashed curves). The size parameter r] = 0.2/ ko and the volume fraction 
fa = 0.25. 
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Figure 7: The angular coordinates a, j3, and the absolute value of the quantity W, plotted versus relative 
size parameter kot] for spheroid orientation angle (p = 90° (blue, dashed curves), 60° (green, solid curves) 
and 30° (red, broken dashed curves). The volume fraction /„ = 0.25 and the eccentricity parameter p = 9. 
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